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We present exact analytical solutions for the zero-energy modes of two-dimensional massless Dirac 
fermions fully confined within a smooth one-dimensional potential V(x) — ~a/ cosh(/3a;), which pro- 
vides a good fit for potential profiles of existing top-gated graphene structures. We show that there 
is a threshold value of the characteristic potential strength aj ft for which the first mode appears, 
in striking contrast to the non-relativistic case. A simple relationship between the characteristic 
strength and the number of modes within the potential is found. An experimental setup is proposed 
for the observation of these modes. The proposed geometry could be utilized in future graphene- 
based devices with high on/off current ratios. 
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I. INTRODUCTION 



Klein proposed that relativistic particles do not experience exponential damping within a barrier like their non- 
relativistic counterparts, and that as the barrier height tends towards infinity, the transmission coefficient approaches 
unity.- This inherent property of relativistic particles makes confinement non-trivial. Carriers within graphene behave 
as two-dimensional (2D) massless Dirac fermions, exhibiting relativistic behavior at sub-light speed^ 3 - owing to their 
linear dispersion, which leads to many optical analogies. 4 >~— Klein tunneling through p-n junction structures in graphene 
has been studied both theoretically^^— and experimentally.— ~— Quasi-bound states were considered in order to study 
resonant tunneling through various sharply terminated barriers^ ^ 13 ! 17 We propose to change the geometry of the 
problem in order to study the propagation of fully confined modes along a smooth electrostatic potential, much like 
photons moving along an optical fiber. 

So far quasi-one-dimensional channels have been achieved within graphene nanoribbonsj 3 "*^— however, controlling 
their transport properties requires precise tailoring of edge termination^ currently unachievable. In this paper we 
claim that truly bound modes can be created within bulk graphene by top gated structures ^LM. such as the one 
shown in Fig. [TJa) . In an ideal graphene sheet at half-filling, the Fermi level is at the Dirac point and the density of 
states for a linear 2D dispersion vanishes. In realistic graphene devices the Fermi level can be set using the back gate. 
This is key to the realization of truly bound modes within a graphene waveguide, as zero-energy modes cannot escape 
into the bulk as there are no states to tunnel into. Moreover the electrostatic confinement isolates carriers from the 
sample edges, which are considered as a major source of intervalley scattering in clean graphene^ 



FIG. 1. (a) A schematic diagram of a Gedankenexperiment for the observation of localized modes in graphene waveguides, 
created by the top gate (Vtg)- The Fermi level is set using the back gate (Vbg) to be at the Dirac point (ep = 0). (b) The 
electrostatic potential created by the applied top gate voltage. The plane shows the Fermi level position at ep = 0. 

In this paper we obtain an exact analytical solution for bound modes within a smooth electrostatic potential in 
pristine graphene at half-filling, count the number of modes and calculate the conductance of the channel. The 
conductance carried by each of these modes is comparable to the minimal conductivity of a realistic disordered 
graphene system^i 3 ^— For the considered model potential we show that there is a threshold potential characteristic 
strength (the product of the potential strength with its width), for which bound modes appear. Whereas a symmetric 
quantum well always contains a bound mode for non-relativistic particles, we show that it is not the case for charge 
carriers in graphene. 



The Hamiltonian of graphene for a two-component Dirac wavefunction in the presence of a one-dimensional potential 
U(x) is 



where <r x>v are the Pauli spin matrices, p x — —iTij^ and p y = —ih-^j are the momentum operators in the x and y 

directions respectively and vf ~ 1 X 10 6 m/s is the Fermi velocity in graphene. In what follows we will consider smooth 
confining potentials, which do not mix the two non-equivalent valleys. All our results herein can be easily reproduced 
for the other valley. When Eq. ([1} is applied to a two-component Dirac wavefunction of the form: 




II. FULLY CONFINED MODES IN A MODEL POTENTIAL 



H = V F (<T x fi x + (JyPy) + U(X), 



(1) 
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where *&a(x) and ^six) are the wavefunctions associated with the A and B sublattices of graphene respectively and 
the free motion in the y-direction is characterized by the wavevector q y measured with respect to the Dirac point, the 
following coupled first-order differential equations are obtained: 

(V(x) - e) 9 A (x) - * + Qy) *b(x) = 0, (2) 



- i (J^ - q y ^j *a{x) + (V(x) - e) * B (x) = 0. (3) 

Here V(x) = U{x)/Tm$ and energy e is measured in units of hvp. 

For the treatment of confined modes within a symmetric electron waveguide, V{x) = V(—x), it is convenient to 
consider symmetric and anti-symmetric modes. One can see from Eqs. (|2]|3|) that ^a(x) and ^b(x) are neither even 
nor odd, so we transform to symmetrized functions: 

*i = V A (x) - iV B (x), * 2 = %a(x) + i^ B {x). 
The wavefunctions and <I>2 satisfy the following system of coupled first-order differential equations: 

d^? 

[V(x)-(e-q y )}^ 1 -—^=0, (4) 



^ + [V(x)-(e + q y )}* 2 = 0. (5) 

It is clear from Eqs. (0][S]) that V^i and ^2 have opposite parity. 

For an ideal graphene sheet at half-filling the conductivity is expected to vanish due to the vanishing density of 
states. When the Fermi energy is at the Dirac point (e = 0) there are no charge carriers within the system, so graphene 
is a perfect insulator. However all available experiments demonstrate non- vanishing minimal conductivity 2 -' 3 -' 3 ^ of the 
order of e 2 /h which is thought to be due to disorder within the system^i - — or finite-size effectsj 34 i 35 In order to study 
confined states within and conductance along an electron waveguide it is necessary to use the back gate to fix the 
Fermi energy (ep) at zero, as shown in Fig. [ljb). Note that Fig. [lja) is just a schematic of the proposed experimental 
geometry and that side contacts may be needed to maintain the 'bulk' Fermi level at zero energy. The conductivity of 
the graphene sheet is a minimum and for a square sample the conductance is of the order of the conductance carried 
by a single mode within a waveguide. Thus the appearance of confined modes within the electron waveguide will 
drastically change the conductance of a graphene flake. Indeed each mode will contribute 4e 2 /h, taking into account 
valley and spin degeneracy. For device applications the sample should be designed in such a way that the contribution 
to the conductance from the confined modes is most prominent. In the ideal case, the conductivity of the channel 
would be the only contribution to that of the graphene sheet. When £p 7^ 0, the conductivity will be dominated by 
the 2D Fermi sea of electrons throughout the graphene sheet. Henceforth, we shall consider the modes for e = 0. 

We shall consider truly smooth potentials, allowing us to avoid the statement of "sharp but smooth" potentials, 
which is commonly used to neglect intervalley mixing for the tunneling problem^— Furthermore, we are interested 
in potentials that vanish at infinity and have at least two fitting parameters, characterizing their width and strength, 
in order to fit experimental potential profiles*^— Let us consider the following class of potentials, which satisfy the 
aforementioned requirements: 

V ^ ~ ~cosh A (/3a;)' ^ 

where a, /3 and A are positive parameters. The negative sign in Eq. ([6]) reflects a potential well for electrons, and 
similar results can easily be obtained for holes by changing the sign of V(x). Notably A = 2 is the familiar case of 
the Poschl- Teller potential, which has an analytic solution for the non-relativistic case^ It is shown in the Appendix 
that the model potential (|6]) with A = 1 provides an excellent fit to a graphene top-gate structure. 

Eliminating ^2 (^1) reduces the system of Eqs. (|4][5]) to a single second order differential equation for ^1 (^2), 
which for the potential given by Eq. ([6]) and e — becomes 

d 2v I / i2 , % wcosh~ A (z) d*i 2 / 9 , i\, N A 9 \ T 

+ Atanh (z) ^' l ' 2 + ^ 2 cosh- 2A (z - A J * 1)2 = 0, 7 

dz 2 wcosh (z)±A dz \ J 
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where we use the dimensionless variables z — (3x, A = q y /(3 and ui = a//3. For A = 1, the change of variable 
£ = tanh(z) allows Eq. ([7]) to be reduced to a set of hypergeometric equations yielding the following non-normalized 
bound solutions for A > 0: 

2' 

4 1 (Si 



#i, 2 = ±(i + p (i-0* 2F1 [p 



ui\2p- 



+ (-!)" {1 + S) q {I- if 2F1 p + g-w, p + q + LJ;2p + 



2 

1 1 



where in order to terminate the hypergeometric series it is necessary to satisfy p = ^ (to — n) + i, q = ^ (uj — n) — j 
and A = uj — n — i, where n is a positive integer. Though we have assumed that uj is positive, one can see that 
the structure of the solutions in Eq. (j8]) remains unchanged with the change of sign of uj, reflecting electron- hole 
symmetry. In order to avoid a singularity at £ = ±1 we require that both p > and q > and obtain the condition 
that ui — n > i. It should be noted that this puts an upper limit on n, the order of termination of the hypergeometric 
series. Notably the first mode occurs at n — 0, thus there is a lower threshold of u> > h for which bound modes appear. 
Hence within graphene, quantum wells are very different to the non-relativistic case; bound states are not present for 
any symmetric potential, they are only present for significantly strong or wide potentials, such that u> = a/ (3 > i. 

Let us consider the first mode (n = 0) in Eq. ([8} which appears within the electronic waveguide with increasing uj. 
In this case the hypergeometric function is unity, and the normalized wavefunctions are: 



where is given by: 



*i, 2 = Ai 2 



(1 + (1 - ± (1 + T+3 (1 - f " 3 



P (2u - l)T {u)T (a + \) 
[2r 2 (uj + i) ± (2w - 1) T 2 (uj)] 



(9) 



(10) 



where T(z) is the Gamma function. As expected, the two functions given by Eq. ([9]) are of different parity, thus unlike 
the non-relativistic case there is an odd function corresponding to the first confined mode. This leads to a threshold 
in the characteristic potential strength uj at which the first confined mode appears; much like in the conventional 
quantum well, where the first odd state appears only for a sufficiently deep or wide potential well. In Fig. [5] we 
present <J/i, ^2 and the corresponding electron density profiles for the first and second bound modes for the case of 
uj = 2. The shape of the confinement potential is shown for guidance within the same figure. The charge density 
profile for these modes differs drastically from the non-relativistic case. The first mode (n = 0) has a dip in the 
middle of the potential well, whereas the second mode (n = 1) has a maximum. This is a consequence of the complex 
two-component structure of the wavefunctions. 




FIG. 2. The wavefunctions (solid line) and ^2 (dashed line) are shown for ui = 2 for (a) the first (n = 0) mode and (b) the 
second (n = 1) mode. A potential profile is provided as a guide for the eye (dotted line). The insets show the electron density 
profile for the corresponding modes. 



For negative values of A , ^1 and ^> 2 switch parity such that $1 2(-A) = \E , 2 ,i(A). This means backscattering 
within a channel requires a change of parity of the wavefunctions in the cc-direction. Notably, when another non- 
equivalent Dirac valley is considered, one finds that there are modes of the same parity propagating in the opposite 



5 



direction. However intervalley scattering requires a very short-range potential or proximity to the sample edges. Thus 
for smooth scattering potentials backscattering should be strongly suppressed. Such suppression should result in an 
increase in the mean free path of the channel compared to that of graphene. This is similar to the suppression of 
backscattering in carbon nanotubes,— where the ballistic regime is believed to persist up to room temperature with a 
mean free path exceeding one micrometer j 38 ' 39 In some sense the considered waveguide can be thought of as a carbon 
nanotube-like structure with parameters controlled by the top gate. 

Notably, our results are also applicable for the case of a one-dimensional massive particle which is confined in the 
same potential. This can be achieved by the substitution \q y \ — > mvp/h into Eqs. (jU[5]), where the gap is given by 
2mvp. Hence in the massless 2D case the momentum along the waveguide plays the same role as the gap in the 
massive one-dimensional case. Therefore, in a massive one-dimensional Dirac system (such as a narrow gap carbon 
nanotube) there exists a bound state in the middle of the gap for certain values of the characteristic strength of the 
potential. 

The number of modes N u at a fixed value of ui is the integer part of w + | herein denoted N u = \u + • The 
conductance of an ideal one-dimensional channel characterized by U) is found using the Landauer formula to be 
G w = AN^e 2 /h. By modulating the parameters of the potential, one can increase the conductance of the channel 
from zero in jumps of 4e 2 /h. The appearance of the first and further confined modes within the conducting channel 
modifies both the strength and the profile of the potential. This nonlinear screening effec t 14 i 15 is neglected in the 
above expression for G u and shall be a subject of future investigation. 

Exact solutions for confined modes can also be found for the Poschl- Teller potential, V(x) = —a/ cosh 2 (/3ir), which 
corresponds to A = 2 in Eq. (J7J. Wavefunctions can be expressed via Heun polynomials in variable £ = tanh(/3x). 



III. DISCUSSION AND CONCLUSIONS 

All the results obtained in this paper have been for a specific potential. However, general conclusions can be drawn 
from these results for any symmetric potential. Namely, the product of the potential strength and its width dictates 
the number of confined modes within the channel.™ Moreover, this product has a threshold value for which the first 
mode appears. The width of the potential is defined by the geometry of the top gate structure, and the strength of the 
potential is defined by the voltage applied to the top gate. The mean free path of electrons within graphene is of the 
order of lOOnm and sub-lOOnm width gates have been reported in the literature^ 9 - - — making quantum effects relevant. 
The number of modes within such top-gated structures is governed by the strength of the potential, with new modes 
appearing with increasing potential strength. In a top gate structure modeled by our potential with a width at half 
maximum of 50nm, the first bound mode should appear for a potential strength of approximately 17mcV and further 
modes should appear with increasing potential strength in steps of 34meV, which corresponds to 395 K. Therefore a 
noticeable change in conductivity should be observed in realistic structures even at room temperature. This is similar 
to the quantum Hall effect which is observed in graphene at room temperature^ A change of geometry, from normal 
transmission to propagation along a potential, allows graphene to be used as a switching device. 

In summary, we show that contrary to the widespread belief, truly confined (non-leaky) modes are possible in 
graphene in a smooth electrostatic potential vanishing at infinity. Full confinement is possible for zero-energy modes 
due to the vanishing density of states at the charge neutrality point. We present exact analytical solutions for fully 
confined zero-energy modes in the potential V(x) = — a/ cosh(/3x), which provides a good fit (see Appendix) to 
experimental potential profiles in existing top gate structures^ - — Within such a potential there is a threshold value 
of u) = a/ (3 for which bound modes first appear, which is different to conventional non-relativistic systems. We found 
a simple relation between the number of confined modes and the characteristic potential strength u. The threshold 
potential strength enables on/off behavior within the graphene waveguide, and suggests future device applications. 
The existence of bound modes within smooth potentials in graphene may provide an additional argument in favor of 
the mechanism for minimal conductivity, where charge puddles lead to a percolation network of conducting channels. 33 

There are experimental challenges which need to be resolved in order to observe confined modes in graphene 
waveguides. These include creating narrow gates and thin dielectric layers as well as optimizing the geometry of 
the sample to reduce the background conductance. Our work also poses further theoretical problems, including the 
study of non-linear screening, many-body effects, parity changing backscattering and inter-valley scattering within 
the channel. 

The study of quasi-one-dimensional channels within conventional semiconductor systems has lead to many inter- 
esting effects. Many problems are still outstanding, including the 0.7 anomaly in the ballistic conductance which 
is a subject of extensive experimental and theoretical study^ 2 - We envisage that the ability to produce quasi-one- 
dimensional channels within graphene will reveal new and non-trivial physics. 
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To illustrate the relevance of our model potential to realistic top-gate structures we provide a simple calculation of 
the potential distribution in the graphene plane for a simplified top-gate structure. Fig. [3] shows the model that we 
use for our estimate. We consider a wire of radius ro, separated by distance h from a metallic substrate (i.e. doped 
Si) and calculate the potential profile in the graphene plane separated from the same substrate by distance d. We 
are interested in the case when the Fermi level in graphene is at zero energy (the Dirac point), hence we assume the 
absence of free carriers in graphene. We also assume the absence of any dielectric layers, but the problem can be 
generalized in such an instance. 



FIG. 3. The simplified geometry used to obtain the model potential. The image charge is shown for convenience. 
One can easily show that the potential energy for an electron in the graphene plane is given by 



where (f>o = <f>a/ln((2h — ro)/Vo), (j)o is the voltage applied between the top electrode and metallic substrate and e is 
the absolute value of the electron charge. One can see that this potential behaves as: 



Appendix: Potential due to a wire above a graphene sheet 





(A.l) 



U{x) 



U + e</> 



2hd 



2 



x < (h-d); 



(h 2 - d 2 ) 2 



x' 



U{x) 




x > (h + d). 



The depth of the potential well is given by 
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and the half width at half maximum (HWHM) is given by xq = y/h 2 — d 2 . In Fig. [4] we show a comparison between 
the potential given by Eq. (jA.ll) and the potential considered in our paper with the same HWHM and potential 
strength. Clearly the potential given by — Uq/ cosh (f3x), provides a significantly better approximation than that of 
the square well potential. 




FIG. 4. A comparison between the potential created by a wire suspended above the graphene plane (solid line) and the 
— Uo/ cosh {fix) potential with the same half width at half maximum (dashed line). 
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